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ELLIPSOIDAL BGK MODEL NEAR A GLOBAL MAXWELLIAN 


SEOK-BAE YUN 

Abstract. The BGK model has been widely used in place of the Boltzmann equation 
because of the qualitatively satisfactory results it provides at relatively low computational 
cost. There is, however, a major drawback to the BGK model: The hydrodynamic 
limit at the Navier-Stokes level is not correct. One evidence is that the Prandtl number 
computed using the BGK model does not agree with what is derived from the Boltzmann 
equation. To overcome this problem, Holway m introduced the ellipsoidal BGK model 
where the local Maxwellian is replaced by a non-isotropic Gaussian. In this paper, we 
prove the existence of classical solutions of the ES-BGK model when the initial data is a 
small perturbation of the global Maxwellian. The key observation is that the degeneracy 
of the ellipsoidal BGK model is comparable to that of the original BGK model or the 
Boltzmann equation in the range —1/2 < v < 1. 
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1. Introduction 

The dynamics of a non-ionized monatomic rarefied gas is governed by the Boltzmann 
equation. But the complex structure of the collision operator has long been a major obstacle 
for theoretical and computational investigation of the Boltzmann equation. To overcome 
this difficulty, Bhatnagar et al.JB], and independently Walender )40j . introduced a model 
equation called the BGK model, where the collision operator is replaced by a relaxation 
operator. Since then, it has been widely used in place of the Boltzmann equation for various 
computational experiments, since this model provides very satisfactory results at relatively 
low computational cost compared to the Boltzmann equation. But the BGK model has a 
major drawback. Hydrodynamic limit at the Navier Stokes level is not satisfactory in that 

Key words and phrases. BGK model, Ellipsoidal BGK model, Boltzmann equation, Kinetic theory of 
gases, Nonlinear energy method. 
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the Prandtl number - defined as the ratio between the viscovity and the thermal conductivity 
- computed using the BGK model is incorrect: The Prandlt number for the Navier-Stokes 
equation is around 0.7, but the computation using the BGK model yields 1. To resolve this 
problem, Holway suggested a variant of the BGK model, called the ellipsoidal BGK model 
(ES-BGK model) [21]: 

d t F + v-S7 x F = A v (M v (F)-F), 

F(x,v, 0) = F 0 ( x,v). 

F(x , v, t) denotes the velocity distribution function representing the number density on the 
phase space point (x,v) in x R® at time t £ R+. A u is the collision frequency whose 
explicit form will be given later. The non-isotropic Gaussian A i v (F) in the r.h.s of (11.11) is 
defined as follows: First, we define the macroscopic density p , bulk velocity U , temperature 


T and the stress tensor 0 by 


p(x,t) = 

/ F(x,v,t)dv, 

J R3 

/* 

p(x,t)U(x,t) = 

/ F(x,v,t)vdv , 

J R 3 

/* 

3 p(x,t)T(x,t) = 

/ F(x,v,t)\v — U(x,t)\ 2 dv, 

Jr 3 

/* 

p(x,t)Q(x,t) = 

/ F(x, v, t)(v — U) <S> (v — U)dv, 

J R3 


and introduce the temperature tensor 71/ as a linear combination of T and 0: 
/ (1 — v)T + vQn 

% = vQ 2 i (1 — v)T + vQ 2 2 

\ ^©31 £''©32 

= (1 — u)TId + £20. 

The non-isotropic Gaussian is now defined as follows: 

We note that the temperature is recovered as the trace of %: 

3 T = ©n + ©22 + ©33 = trQ = trT v . 


1/013 

££©23 

(1 - v)T + vQ 33 


The collision frequency A v takes the following explicit form: 


A v 


pT 

1 - £/’ 


-2 <Z/< 


1. 


The free parameter v is introduced to derive the correct Prandtl number. The restriction 
on the range of v is imposed to guarantee that the temperature tensor Tv remains positive 
definite. (See [2]). Prandtl number computed via the Chapman-Enskog expansion using the 
ES-BGK model is given by Pr = 1/(1 — v) (See [2] [TUI ECO EH)- The two most important 
cases in the range —1/2 < v < 1 are v = 0 and v = (Pr— 1)/Pr « —3/7: When v = 0, (11.11) 
reduces to the classical BGK model, whereas v = (Pr — 1 )/Pr corresponds to the ES-BGK 
model with the correct Prandtl number. 

The relaxation operator of the ES-BGK model satisfies the following cancelation property 
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mm- 


I (M„(F)-F) 

Jr 3 



dv = 0, 


which implies the conservation of mass, momentum and energy: 



F(x, v, t)\v\ 2 dxdv = I F 0 (x,v)\v\ 2 dxdv. 



Entropy dissipation property was established recently in [2] : 


7 /* 

— / F(t) \ogF{t)dxdv < 0. 

™ Jt 3 xR 3 


It is important to note that, as in the case of the original BGK model or the Boltzmann 
equation, the only possible equilibrium state for EU is the local Maxwellian: 



To see this, let’s assume that A!„(/) = /• We then recall the definition of 0 and T v to see 
that 



j M v {F)(v){v - U) ® (v - U)dv = p%. 

Therefore, upon multiplying (v — U) ® (v — U ) to both sides of M V (F) = F and integrating 
with respect to v, we have 

p% = p&. 

In view of the definition of T v , this leads to 

(1 - v)TId + vQ = 0. 

Thus, 0 = TId , and we see, from the definition of that T v = TId for 3x3 identity 
matrix Id. This gives 



= - P 2 T 

(2ttT) 3 /2 

= M(F). 


P \v-uf 


That is, Ad v [F) reduces to the usual local Maxwellian M(F). 
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In this paper, we study the existence of classical solutions of (O) and their asymptotic 
behavior when the initial data is a small perturbation of the normalized global Maxwellian: 

(L3) 

We define the perturbation / around /r by the relation: Fix , v , t) = ^l + y/JIf(x, v, t) and, ac¬ 
cordingly, -£ 0 ( 2 ;, v) = n+y/Jifo(x, v). Then, after linearization around the global Maxwellian, 
the ES-BGK model takes the following form (See Section 2 for precise definition of each 
term): 

&tf + v-V x f = L v f + T(f), 
f{x,v,0) = f 0 {x,v). 

where L„ denotes the linearized relaxation operator and T(/) is the nonlinear part. In section 
2, we verify that L v can be represented as a ^-perturbation of the linearized relaxation 
operator of the original BGK model: 

(1-5) L v f=(P Q f-f) + vP 1 f + vP 2 f. 

Here, P 0 denotes the macroscopic projection operartor on the the linear space generated by 
{y/JI,Vy/jI, | Pi and P 2 are operators related to the burnett functions, which play 
a crucial role in the hydrodynamic limit of the Boltzmann equation at the Navier-Stokes 
level. (See |3]). In general, the coercivity estimate of the linearized collision or relaxation 
operators for spatially inhomogeneous collisional kinetic equations are degenerate, and the 
major difficulty in obtaining the global existence in the perturbative regime lies in removing 
the degeneracy to recover the full coercivity T8j HU [20]. When the spatial variable lies 
in T 3 , the usual recipe is the use of the Poincare inequality together with a system of 
macroscopic equations and the conservation laws (See, for example, [HE ED])- In the whole 
space, where the Poincare inequality is not available, additional consideration has to be 
made to compensate the still lingering degeneracy pea hu m mi ng eg]. Therefore, it is 
very important to capture the degenerate coercivity estimate of the linearized relaxation 
operator first. In our case, it is not clear whether the presence of the additional terms Pi 
and P 2 make the linearized relxation operator more degenerate or not. In Theorem 12.81 we 
show that, for —1/2 < v < 1, the degenerate coercive estimate of L„f is comparable to that 
of Lof = (Po — /)/, for which the usual energy method is well-established (See Theorem 

EE}: 


{L v f,f) L % < -a||{I-P 0 }/||i;, (-1/2 <«/<!), 


for some constant C v > 0. This indicates that the dissipative property of the linearized 
relaxation operator for the ES-BGK model is essentially same as that of the BGK model or 
Boltzmann equation. 

On the other hand, since the ES-BGK model is obtained by replacing the temperature 
function T by the temperature tensor 7/ in the classical BGK model, additional difficulties 
related to 7/, which was not observed in the classical BGK model arise. First, in each step 
of the iteration scheme designed to obtain the local in time existence of the solution, we need 
to check that the temperature tensor remains strictly positive definite, which is established 
in Proposition lh.il as: 


C~ 


T{x,t) 


Id<T v 1 (x,t) < 


C~ 


T(x,t) 


Id. 
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where C v \ = min{l — v, 1 + 2z/} and C v \ = max{l — v, 1 + 2u}. This also shows why 
the restriction of the range of the free parameter in the interval (—1/2,1) is crucial: It is 
only in this range that the temperature tensor is comparable to T, and, therefore, the non¬ 
isotropic Gaussian is comparable to the local Maxwellian. Secondly, due to the presence of 
the free parameter v in the definition of the temperature field 7/, it is a priori not clear 
whether the nonlinear perturbation T(/) can be estimated uniformly with respect to v near 
v = 0 because the the inverse of the temperature tensor 7/T 1 may have problematic terms 
involving 1 jv. Such a singularity at v = 0 is undesirable considering that the case v = 0 
corresponds to the classical BGK model. The above equivalence estimate guarantees that 
such singularity never shows up when —1/2 < v < 1. 

The mathematical theory for the BGK model has a rather short history. The first rigor¬ 
ous existence result can be traced back to Ukai m , where he considered stationary problem 
for 1 dimensional BGK model in a periodic bounded domain. Perthame established the exis¬ 
tence of weak solutions of the BGK model with constant collision frequency in m assuming 
only the finite mass, momentum, energy and entropy. See also [7]. The uniqueness was con¬ 
sidered in a more stringent functional space involving the pointwise decay in velocity 128] . 
Mischler considered similar problems in the whole space in |25l . Extension to L p was carried 
out in [42]. Issautier established regularity estimates for the BGK model and proved the 
convergence of a Monte-Carlo type scheme to the regular distribution function in [22]. The 
convergence property of a semi-Lagrangian scheme for the BGK model was studied in [29| . 
In near Maxwellian regime, Bellouquid [5] obtained the global well posedness in the whole 
space using Ukai’s spectral analysis argument [36]. In the periodic case, Chan employed the 
energy method developed by Liu et al. [24] to establish the global in time classical solution 
near global Maxwellians m- The convergence rate to the equilibrium was not known in 
this work, which was derived in [41] . For fluid dynamic limit of the BGK model, see [301131] . 
The ES-BGK model has attracted only limited attention until very recently since it was not 
clear whether the entropy dissipation property holds for this model. It was proved in the 
affirmative, at least at the formal level, in [2], which revived the interest on this model. 
To our knowledge, no existence result has been established for the ellipsoidal BGK. For 
numerical test for the ES-BGK model, we refer to p] muni Esiiii]. For general review of 
the mathematical and physical theory of the Boltzmann equation and the BGK model, see 

pnunnsiisiEsusg. 

Before proceeding further, we define some notations. 

• When there is no risk of confusion, we use generic constants C. Their value may 
change from line to line but does not depend on important parameters. 

• We define the index set i < j by 

Clij = di 2 + 023 + O 31 . 

i<j 

• ei (i = 1,2, 3) denote the standard coordinate unit vectors in M 3 . 

• 0” denotes n-dimensional zero vector. 

• /(to, n; a, b) denotes a (m + n) x (to + n) diagonal matrix whose first to diagonal 
elements are o and following n diagonal elements are b. 
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• (-,-) i 2 and (-,-) i 2 denote the standard L 2 inner product on R^ and x 
respectively: 


(L9)lz = f f(v)g(v)dv, 

Jr 3 

(/:ff)L 3 „=/ f{x,v)g(x,v)dxdv. 

' J T 3 xR 3 


/T 3 xR 3 

• || • 11^2 and || • ||i 2 denote the standard L 2 inner norms on R^ and x R^ respec¬ 
tively: 


\h 


2 = ([ \f(v)\ 2 dvV, 
v Jr 3 7 

= ( [ \f{x,v)\ 2 dxdvY. 
v ' JT 3 xl 3 7 


l£; 


• We employ the following notations for the multi-indices and differential operators: 

a = [a 0 ,011,0:2,0:3], /3 = [^i,/3 2 ,/3 3 ], 

and 


a? = d?°d%d%d%d^d%d%. 


For simplicity, when only the spatial derivatives are involved, we write 9“ = 9“* 9“ 2 2 9“ 3 . 

1.1. Main results. We now state our main result. We first define the high order energy 
functional £(f(t)): 

£(m) = \ e \\ d pm\\i lv + E tmmwb^ds. 

\a\ + \/3\<N ' \a\+\0\<N J ° 

Theorem 1.1. Let — 1/2 < v < 1 and N > 4. Let Lo = // + i/Ti/o > 0 and suppose 
fo satisfies i2.5\) . Then there exist positive constants 5 V and C = C{N,v), such that if 
£(fo) < then there exists a unique global solution f to o such that 

(1) The distribution function is non-negative for all t > 0: 

F = n + yffif > 0, 

and satisfies the conservation laws m - 

(2) The high order energy functional £ (/(i)) is uniformly bounded: 

£(m)<c£(f 0 ). 

(3) The distribution function converges to the global equilibrium exponentially fast: 

E \\d%m\\ Llv <Ce- c,t 

M + |/3|<iV 

for some constant C and C'. 

(4) If f denotes another solution corresponding to initial date fo satisfying the same 
assumptions, then we have the following uniform L 2 -stability estimate: 

ll/(i)-/WIUl,„<C|l/o-/o|Ui ! „. 
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This paper is organized as follows: In section 2, we consider the derivation of the linearized 
ES-BGK equation and the main result is stated. We also derive the coercive estimate and 
determine the kernel of L„. In section 3, various estimates on the macroscopic field are 
established and, based on this, the local in time existence is obtained. In section 4, the 
nonlinear energy estimate is derived, which readily leads to the global existence and the 
asymptotic behavior. 


2. Linearization 

In this section, we consider the linearzation of the ES-BGK model around the global 
Maxwellian m- For some technical reason, we define G u as follows: 


„ 1-v [3pT + p|C/| 2 

Gv -~r\ 2 


Id + u[ pe + f 9 U )-^Id. 


Due to the symmetry of G , we can view G as an element in 


{G n, G 22 , G 33 , G 12 , G 23 , G 31 } . 


We also define J v to be the Jacobian matrix for the change of variable (p, U, %) —► (p, pU , G): 


_ d(p, pU, G v ) 

d{p,u,%) ■ 

Lemma 2.1. (1) J v is given by 


( 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

\ 

Gi 

9 

0 

0 

0 

0 

0 

0 

0 

0 


u 2 

0 

p 

0 

0 

0 

0 

0 

0 

0 


u 3 

0 

0 

p 

0 

0 

0 

0 

0 

0 



±f^pU x 

^pu 2 

^pU 3 

\p 

0 

0 

0 

0 

0 


A v 

^22 

^ PUi 

^pU 2 

^pU 3 

0 

\p 

0 

0 

0 

0 


^33 

^PU 1 

^ P u 2 

l -^pU 3 

0 

0 

\p 

0 

0 

0 


/1 12 

IPU'2 

1PU 1 

0 

0 

0 

0 

\p 

0 

0 


^23 

0 

%pU3 

1 PU 2 

0 

0 

0 

0 

\p 

0 


V ^31 

iPU'i 

0 

1 PU 1 

0 

0 

0 

0 

0 

\p 

/ 


where is defined as 
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(2) J v x is given by 











( 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

\ 


_(7i 

1 

0 

0 

0 

0 

0 

0 

0 

0 



_U 2 

0 

3_ 

0 

0 

0 

0 

0 

0 

0 




0 

P 

0 

_1 

0 

0 

0 

0 

0 

0 



An 

2 (l+ 2 v) U-! 

3 p 

2(1-!/) (7 2 
/ 3 s P 

2(1-!/) (73 

, 3 , p 

2 

P 

0 

0 

0 

0 

0 



A 22 

2(1-!/) Ui 

2(l+2i/) f/ 2 

2(1-1/) u 3 

0 

2 

0 

0 

0 

0 



A33 

2(1-!/) t/l 

3 p 

2(1-!/) C/ 2 

3 P 

2 (l+ 2 v) Us 

3 P 

0 

0 

2 

P 

0 

0 

0 



A12 


-i/a 

0 

0 

0 

0 

2 

0 

0 



A23 

0 

-i/a 

-i/a 

0 

0 

0 

0 

2 

0 



\ A31 

-i/fia 

P 

0 

-A 

p 

0 

0 

0 

0 

0 

2 

P 

/ 

where Af- is 

defined as 












- + 

(l^)|E/| 2 

+ 3ot/f 

-il 

l. 





Aij — ( Tij + UiUj). 

(3) When F = p, J v and J~ x reduce to the following simpler form: 

Jv \ F=ft = 1(4,6; 1 , 1 / 2 ) and J~ l | F=/i = 1(4, 6 ; 1 , 2 ). 

For i/ie definition of I(m,ri',a,b), see the notation at the end of the introduction. 

Proof. The proof is straightforward but very tedious. We omit the proof. □ 


Lemma 2 . 2 . We have 

(1) Derivatives for det7D.‘ (1 <i,j < 3, i ^ j) 


ddetTv 

dp 

d det 70 

arili 


= 0 , 

f=h 

= 1 , 

F=v 


d det 71/ 
dUi 

d det 71- 

^T 


= 0 , 

= 0 . 


( 2 ) Derivatives for M v : (1 < i,j 


dM v 

dp 

dM v 

dTu 


F=n 

F=fM 


— 3, i 7 ^ j) 


p(v), 


dMv 


dUi 


F=n 



= Vip(v), 



ViVj)u(v). 


Proof. (1) A straightforward calculation leads to the following explicit form of the determi¬ 
nant of Tv- 


(2.1) det70 — 7n 722*733 — T^fTii — T 31 T 22 — T^Tn- 

Then (1) follows from explicit calculations using 

dJf L = n d JfL = n S = / 

a P ’at/ ’ ar ffe 1 


and 


i ( i = e : j = k ) 

0 ( otherwise ) ’ 
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(2) We only consider 4^-. Other terms can be obtained similarly. We first observe that 


dM, 


dTij 

When i = j = 1, we have 
dM „ 


1 1 <9detT 1, . i 

2 |detT| ar, + 2^ _ 


ar 

8 %, 


T~ 1 (y - U) 


M v , 


arn 


F=/x 


= < —- + -v 


2 2 


(1 

0 

0 ^ 

1 

0 

0 

0 

v 

V o 

0 

0 ) 

1 j 


v\ — l 


P- 


can be obtained in the same manner. In the case i ^ j, we observe that 

dM , 


Similarly, we have 


dTu 


dM 



0 + iu T 

f 0 

1 

0 ^ 

1 

= 


0 

0 

V 

a. 

II 

fc, 

2 

^ 0 

0 

oj 

1 


dT 23 


= fl^/T. 


, dM 

= v 2 v 3 p, and —— 
F—n dT 3 1 


F=/x 


= V 3 V\H- 


□ 


Now, we are ready to prove the main theorem of this section, which basically says that 
the linearized relaxation operator is composed of ^-perturbation of the projection on the 
macroscopic kernel and nonlinear terms. 

Theorem 2.3. Let F = fj, + y/JIf . Then the ellipsoidal Gaussian M^(F) can be expanded 
around fj, as follows: 

M„(F) = p + (P u f)y/JI+ J2 ( l\ D t P e, P eU s ,G s) M(e)} (l^efde) 

l<i,j < 3 V-'o J 

Here, P„ is given by a v-perturbation of the usual macroscopic projection Pq : 

Pvf = Pof + v{P\f + P 2 f), 

where 


3 vf - M 


Pof = fy/jidvj v^I+ fvy/jldvj ■ Vy/Ji 

«/=E(//- 

P 2 f = 5Z (^y f’UiVjy/Jidvj ViVjy/ii. 


f^L^tfidv 


l<-3 

V6 




3 V 2 


■\/m, 


i<j 

and M v {&) denotes 


M V {B) = 


Pe 


yJdet(2nTg) 


ex P ~U e ) T g 1 (v - U e ) 
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where the transitional macroscopic fields pg, Ug, Gg and Tg are defined by 
pg = 0p + (1 — 9 ), pgUg = dpU , and Gg = 9G , 

and 

/ (1 - v)Tg + vQgu 

Tg = I 6i2i 

\ ^0031 

Proof. We define g{9) as 

9{0) = M(9(p,pU,G) + (l-9)(l,0 3 ,0 6 )) 

= A i (pg, pgUg,Gg) . 

Note that g{9) represents the transition from the global Maxwellian p(v) to the ellipsoidal 
Gaussian A4^(F). Then we have from the Taylor’s theorem 

(2.2) g( 1) = g( 0) + g'( 0) + g” (9)( 1 - 9) 2 d9. 

Jo 

The first term in the right hand side is the global Maxwellian: g(0) = p. We now consider 
the second and the third terms: 


I/06U2 ^0013 

(1 — v)Tg + 1/0022 ^0023 

^0032 (1 — v)Tg + 1/0033 


(i) g'(0): We observe from Lemma T2.21 that 

n \ a /n\ (dM v dM v dMv\ 


( 


= l,f, 


F=g. 


V7 — 1 Vo — 1 vi — 1 


-,V 1 V2,V 2 V 3 ,V 3 Vi 


>)■ 


2 ’ 2 ’ 2 

Then, using the identities in Lemma HZTT1 and Lemma T2.21 g'( 0) can be represented as 


s'(0) = —M(9(p, pU,G) + (1 — 0)(1,0,0,0,0 6 )) 


(p-l,pU,G) r - Jg-'D^^M 
{p— 1, P U 7 G) T ■ J~ l 


0=0 


?=o 


X | l,v, 


v\ — l v 2 — l v 2 — l 


-,V 1 V2,V 2 V 3 ,V 3 V 1 P 


J f^fpdv^j p + {^J fv^pdv^j ■ vp 




vf - 1 


p T 2 ^ ) GijViVjp. 

i<j 


Here Jg denotes ^fpJijgTe) an< ^ we usec ^ ^o = 3 ■ 
(ii) g" (9): By an explicit computation, we find 


g"(0) = 


d 2 M 

d9 2 


(9(p ~ 1, pU,G) + (l — 9)(1, 0,0 6 )) 


= (p — 1) pU,G) T \ D 


\ Pe , PoUo , Go ) 


M(9) > (p — 1, pU,G). 


(iii) We claim that 


ff'( 0 ) = Pvfyjp■ 
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Note that it is enough to establish 

3 ' „2 




V6 

3 2 ' 3«?-M 2 


£ f 


3x72 




3 y i - H : 

3x/2 


■VJd- 


We first observe that Gu (i = 1,2,3) can be decomposed as 


Gu = 


1 - v 


' f dv + v I f 

R 3 A 

v \ 2 — 3 


I*- r i 

2 Jr 3 2 


6 


3-uf - M 


fdv 


Viddv, 


so that 
3 


u 


vf - 1 


fGu . 

R 3 1 


V^= 2 E / / 


2=1 

3 


\v \ 2 - 3 (3vf- \v\ 

-:-r V 


2=1 

3 


M - 3 


= 2 E( / 


X; r — 1 

VJddv ) 


+ 2 ‘-E(Lj - k ' 1 ' 1 

2 — 1 

+ A + B. 


R 3 


6 


v? - 1 


VJd 


We compute A as 


A = 




6 


2=1 


vf - 


) Vm 


For B , we observe that 


and 




vf - 1 3t>? - M 2 Id 2 - 3 


6 


^21aH_ = o, 


to derive 


6 


5 = 2 E( / f 3V * J V]2 Vfddv 


2=1 

3 


R 3 


2 = 1 
3 


3^ 2 - |^| 2 , H 2 - 3 \ 




6 6 

VJd 


Vt 
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= E 


i— 1 

2 


3vf — \v\ 2 
f —ttts— V^ dv 


/E 


3 V 2 

3v 2 - 


% v i ~ M : 

3y/2 


Vt L 


2 = 1 


6 


y/pdv^j 


- 3\ 


= E 


2=1 


=0 

3v 2 - \v\ 2 


f 


^Jfidv 


3 y j ~ M : 

3y/2 


J 


/r3 3y/2 

Plugging (i), (ii), (iii) into 12.21 we obtained the desired result. 
We now consider the linearization of the collision frequency. 




□ 


Lemma 2.4. The collision frequency A u can be linearized around the normalized global 
Maxwellian as follows: 


A u = 


1 


1 


where 


Ap — 


J J(f l {T e , 0 3 , l/3p e Id)(l - 0)d0^ • (p - 1, P U, G). 
Proof. We expand A v by the Taylor’s theorem. Then the second term reads 

D (petPeUe , Ge) (p 9 Tg){ 1 - Q)dB^ ■ {p — 1, P U,G). 


Note that 
to see 


D (P9,PeUo,G e ) ~ Jq 1 D(p 0 ,U 0 ,Te) > 


Art — 


f Jg^DipePe.TeMm 1 - 0)d0 j • (p-l,pU,G) 
j Ji 1 {T e ^,l/3pgId){l-e)d9^ -(p-l,pU,G). 


□ 


Instead of writing down D 2 pg pgUg G ^ explicitly, we introduce generic notations which 
considerably simplify the argument. We first observe that 

= ■?i lD (pe,Ue,Te)Ji lD (pe,Ue,Te)M(d)- 
We then invoke Lemma f2. II to conclude the following lemma. 

Lemma 2.5. There exist generic polynomials Pff, Rf^ such that 

(p-i, p u,g) t {rP lpt , ptU „ Gl) M(e)} { P — 1 , P u, g) 


= E 




P^( pg,V-Ug,U e ,r g - \lz) 
RM (pg, det Tg) 


exp ( -i(u- Ug) r T e 1 {v~Ug) ) ( f,ei){f,ej ) 


where P/j(x i,..., x n ) and Rfj(xi, ..., x n ) satisfy the following structural assumptions {Rm) 
• (TLm„ 1) Pij u Is a polynomial such that Pij( 0, 0,..., 0) = 0. 


hj 

(Rm„ 2) Ri^A is a monomial. 
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Lemma 2.6. There exist generic polynomials Pf v , Rf v such that 

{j e -\T e ,o 3 ,i/3 Pe id)} ■ (p-i,pu,G) = —</.*>, 

where Pjf ( x\,...,x n ) and Rfj(xi,...,x n ) satisfy the following structural assumptions 

(n A „). 

• (TLa v 1) Pf" is a polynomial such that Pij( 0, 0,..., 0) = 0. 

• {TLa u 2) Rf" is a monomial. 


Note that P'ff, Rif, Pf 4 and R:f u are defined generically. They may change line after 
line during the argument. But explicit form is not important as long as we keep in mind 
the structural assumptions TLm and TLa v ■ To simplify the notation further, we define Qff 
and Qf v as 


and 


3 _ f 1 P% v {pe,v-Uo,Tf\v) 

o Rff{p g ,de t%e) 


exp ( --(« - U e ) 1 T~ e \v - U e ) ) (1 - 9) 2 d0 


/■ 1 p A , 

Qf»(9)= / -A- 
JO 


P^{pe,U e ,Te,v) 


(1 - 6)d6. 


lo Rf^iPe) 

Then the relaxation operator and the collision frequency can be expressed in a more succinct 
form: 

M„(F)-F={P„f~f)^jI+J2 Q% v </> e i)Ll </, ef) Ll , 

and 

^ E e ^?- 

1 — 1/ 1 — 1/ ” 

We summarize the result in the following proposition. 


Proposition 2.1. The relaxation operator can be linearized around the normalized global 
Maxwellian p as follows 

K (M V {F) - F) = E- (l + E Of" if, e<)) { (Puf ~ f) + E 2E </> *)</> e,->} 

» i,j 


We now substitute the standard perturbation F = p + yfpf into (1.1) and apply propo¬ 
sition [2T] to obtain the perturbed ES-BGK model: 


(2.3) 


dtf + v- V x f = L v f + T(f), 
f(x,v, 0) = f 0 (x,v), 


where fo(x,v) = The linearized relaxation operator L u and the nonlinear perturba¬ 

tion T(/) are defined as follows: 

(2-4) Kf = ^—{Pvf~f), 

1 — V v 7 

and 

r (/) = J^{T,Qt(Lei)}(P,f-f) 

i 

+ E 

l<i,j<3 
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+ l v y, Qi e i)Ll(fi e j)Llif, ek) L l 

1<M<3 

= ri(/,/) + r 2 (/,/) + r 3 (/,/,/). 


The conservation laws in (11.21) now take the following form: 

/ f(x,v,t)^ffi dxdv = / fo{x, v)^ffi dxdv, 

-/t|xR3 J T|xRg 

f(x, v, t)vy/JI dxdv = / fo(x,v)v^/JI dxdv, 
Jfl xR?. 


/TjxR?, 


/ 


tj xe; 


f(x, v, t)\v\ 2 y/JI dxdv = 


/T3xR3 


f 0 (x,v)\v\ 2 ^/jl dxdv. 


Therefore, if initial data shares the same mass, momentum and energy with /r, the conser¬ 
vation laws read 


/T3xR3 


f(x, v, dxdv = 0 , 


( 2 . 5 ) 


TJ X R3 


f{x,v,t)v^ffi dxdv = 0 , 


TJxRJ 


f(x, v, t)\v\ 2 ,/Jl dxdv = 0 . 


2.1. Analysis of the linearized relaxation operator. We now study the dissipative 
mechanism of the linearized operator. We start with the following technical lemma. 

Lemma 2.7. Po, Pi and P 2 satisfies the following properties: 

(1) Po, P 2 and P 2 are orthonormal projections: 

P 2 = Po, Pi = Pi, P 2 2 = P 2 . 

(2) Po, Pi and P 3 are mutually orthogonal in the following sense: 

P 0 P 1 = P 1 P 0 = P 0 P 2 = P 2 P 0 = P 1 P 2 = P -2 P\ = 0 . 

Proof. (1) The first and third identities Pq = Po and Pf = P 2 follow from the fact that 
{^/Ji, Vy/JI, \v\ 2 y/JI} and {v\v 2y /fl, v 2 v 3 ^/Ji, v 3 v\yffi} form orthonormal bases respectively. To 
show P 2 = Pi, we first observe that 

<( 3^ 2 - \v\ 2 )VP, ( 3v i - \v\ 2 )y/fi) L2 = 12 , (1 = 1 , 2 , 3 ) 

((3v 2 - \v\ 2 )^/Jl, (3vj - \v\ 2 )^JI) l I = -6 (i ^ j). 

Using this, we have for cfiv) = (3i > 2 — |d| 2 )/3v / 2 
P lf = Pl{(f, Cl) L 2d + (/, C 2 ) L 2C 2 + (/, C 3 ) L 2C 3 } 

= ^{2(/,Ci) L 2 - (/,C 2 ) L 2 - ( f,C 3 ) L 2}Cl 

+ + 2(/, C 2 )x,2 — (/, C 3 ) L 2 }ci 

+ Ci)i ,2 - (/, c 2 ) L 2 + 2(/, c 3 ) L 2 }ci 


= / 


2 ci — c 2 — c 3 


)l ? Ci + ( / - 


-Cl + 2c 2 - c 3 


3 /Li V’ 3 

= {f^C^Lid + (/,C 2 ) L 2C 2 + (f,C 3 ) L 2C 3 


) L ? c2+ ( f - 


-Cl - c 2 + 2 c 3 \ 

- 3 - 


c 3 


3 
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= Plf- 

In the last line, we used c\ + C 2 + C 3 = 0. 

(2) Straightforward calculations gives 

(a/M, (3 v i - \v\ 2 )^JI) L 2 = (3 vf - \v\ 2 )y/fZ) L 2 = ((|u | 2 - 3 (3v 2 - |i;| 2 ) v / m)lj = 0 , 

and 

(viVjy/Ji, (3 v\ - \v\ 2 )^JI ) L 2 = 0. 

This implies (2). □ 

We now prove the main theorem of this section. Note that that the estimate is uniform 
with respect to v. 

Theorem 2.8. For —\<v<l,we have 

( L v f,f ) L 2 < - min jl, ^ }||(/~ Po)f\\ 2 L 2 . 

Proof. From the definition of L^, we have 

(1 - v){L v f, f) L 2 = ( P v f - /, f) L 2 
( 2 . 6 ) =(Pof-f + v(Pi+P 2 )f,f)Ll 

= -\\(I-P 0 )f\\l l +u((P 1 +P 2 )f,f) L 2. 

We recall from Lemma \‘2. 71 that (Pi + P 2 ) _L Pq , which gives 

((P 1+ P 2 )f,f) Ll = ((Pl+P 2 )(I-P 0 )f,(I-Po)f)Ll 

= \\(Pl + P 2 )(I ~ P 0 )f\\ 2 L 2 v . 

We then observe from Lemma o that Pi + P 2 is a projection operator: 

(Pi + P 2 ) 2 = P 2 + PiP 2 + P 2 Pi + P 2 = Pi + P 2 , 

which leads to 

(2-8) ||(Pi + P 2 )(/ - P 0 )fhl < ||(/ - Po)f\\h- 

Therefore, we have from (12.61) - (12.71) 

(1 - v){L v f, f ) L 2 < - min{(l - v), (1 - M)}|| (/ - P 0 )/||£;- 
Since (1 — u) >0, this completes the proof. □ 

Corollary 2.1. For —1/2 < v < 1, the kernel of the linearized relaxation operator is given 

by 

Ker{L v } = Ker{L 0 } = span{yfjl, \v\ 2 ^/Jl}. 
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3. Estimates on the macroscopic field 


3.1. Estimates on the macroscopic field. To control the nonlinear perturbation T(/) 
in the energy norm, we first need to establish various estimates for macroscopic quantities. 
Throughout this section, C„ > 0 means that C v is strictly positive for all — 1/2 < v < 1 . 


Lemma 3.1. Let Eft) be sufficiently small, then there exists a positive constant C > 0 and 
C v > 0 such that 

( 1 ) \p(x,t) - 1 | < cVW), 

(2) | Ui(x,t)\<Cy/£(fj, (l<i<3) 

(3) \T v u(x,t) — 1| < C v y/sjf), (l<t<3) 

( 4 ) \Tij(x,t)\ < uCy/£{L), {l<i<j< 3) 

Proof. (1) We have from Holder inequality 

\ P {x,t)-l\= [ fy/JIdV < \\f \\ L 2 < y/Ift). 

JR 3 

(2) Using the lower bound estimate of p , Holder inequality and f R3 pvdv = 0, we see that 


m = - 

p 


< 


I fviy/pdv 

R 3 

VW) 


< 


LI 


(3) For the upper bound of Tu 


1 


i -Vm 

V ^_ < Cy/£(fj. 

-VW)~ 


Tu = (1 - v)T + vQa 

+ 2 }+4-/ (p + Vhf)v-dv-uA 

. P Jr 3 J l P Jr 3 J 


< 


< 


3 

3 P 

1-u ( 3+V^5||/||l 2 


3+ [ f\v\ 2 y/Jidv\ + — /1 + [ fvjyffidv\ 
Jr 3 J Pi Jr 3 J 


i + v^ll/IUi „ 


< 

< 

Therefore, 

(3.1) 


i + CMli, 

P 


1 + q#) 

1 -V£(t) 


Tu - 1 < < Cy/£(f). 


i -Vm 

Using the lower bound estimate for p and Ui , we estimate the lower bound similarly as 


Tu = (1 - v)T + vdu 

= 0 + Vhf)\v\ 2 dv - |U| 2 } + - { / (p + VPf)vidv-U?\ 

op {Jr 3 J P {Jr 3 J 

= ^ (3+ f f\v\ 2 ^pdv-\U\ 2 X+ -ll+ f fv*y/jIdv-U?\ 

3p l Jr 3 ) P { Jr 3 ) 

> ^ {3 - \/15||/||lJ i „ - C£(t)j + - {1 - v^ll/IU^ - C£(t)} 
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l-C v \\f\\ Ll ' V -CS(t) 

p 

l-C v ,j£(f)-CE{t) 

1 + V£(t) 

1 

i + Vs{t) 


Hence we have 


(3.2) 


Tu-l> 


-C„VE{t) 

1 + yj £{t) 


> -C v yJW). 


(ED and (13.21) give the desired result for Ta (i = 1,2,3). 
(4) Tij can be estimated similarly as 


\Tn\ < 

< 

< 

< 


fviVjy/Jldv + v\Ui\\Uj 




1-JW) 

v\j£{t) 

1 ~VW) 

vCyJW). 


vC£{t) 


;C£{t) 


□ 


Lemma 3.2. Suppose £(t) is sufficiently small. Then there exists a positive constant C\ a \ > 
0 and C\a\,v > 0 such that 

(1) \d a p(x,t)\ < \/£{£), 

(2) \d a U(x,t)\ <C w £(t), 

(3) \d a Tij{x,t)\ <C WUv £{i). 


Proof. (1) Since d a p = 0, we have 


\d a p\ = 

< 


+ = j 


\d a fhi„ < vm. 

(2) A straightforward computation using U = ^ J fv^/JIdv and the chain rule gives to 

\d a U\<^( E [ \d ai f\\v\y/pdv\ (l+ E I d a2 p\ 

P \ I „ I I \ \ 


|ai[:<JV 


|a 2 |<IV 


Then the use of Holder inequality and the estimate (1) leads to 

M 


\d a u\ < q Q| ^EKi±v^E_ < c H VW)- 


(i -m) 
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(3) Recall T, 3 = {/ f\v\ 2 ^/Jidv — \U\ 2 } + ^ {/ fvfy/JIdv — U 2 }. Therefore, by the same 

argument as in ( 2 ) above, we have 


\d a r tJ \ < 


a 


\a\,iy 


9 2|a| 


E / fi+ E +C\a\,v£(t) 

K \ cn \< N jR3 ) \ | Q2 |<iV J 


< c vw)(i+vm)' 1 c 

< c lal „vm- 


□ 


Lemma 3.3. Let Eft) be sufficiently small. Then, we have positive constants C > 0 and 
C v > 0 independent of 9 such that 

(1) \pg(x,t) -1| < Cy/S®, 

(2) \Ug(x,t)\<Cy/S^j, 

(3) I Teu(x,t) - 1| < C v ^/W), (i = 1,2,3), 
{A)\T ei j{x,t)\<C v y/£$), 

Proof. (1) By Lemma 13. II and the definition of pg , we have 


I Pe - 1| = 9\p - 1| < Qy/Etf) < ^£ff). 

( 2 ) follows directly from ( 1 ), Lemma I5TT1 and the definition of Ug: 

Ug = —pU. 

Pe 

(3) We divide the case into i = j and i ^ j. When i = j, we have from the definition of Gg 
that for i = 1 , 2 ,3: 

1 — v ( Pe(®9ii + 0022 + © 033 ) + Pe\Ug\ 2 \ (PeQeu + pUgi\ pg 


(3.3) 


= 9 


1 — V / p(Q 11 + 022 + © 33 ) + p\uy 


pOu + pU\ 


Summing over i = 1,2, 3, we obtain 
(3.4) 


Pe(Qe 11 + 0022 + 0033) _ 0^(011 + 022 + © 33 ) 


, Pe\Ug\ 2 — 9p\U\ 2 , 3, 

+- 2 - + 2^ P6 ~ P °'' 


We substitute (13.41) back to (13.31) to get 


( P9&9U + peUg t \ _ _ 




(3.5) 


1 - v 


P(011 + 022 + 033) + |t^| 


2 \ "1 


1 - V 


{Pe ~ dp) 


9 


\ — V ( p(011 + 022 + 033) + \U\‘ 


+ 


3 V 
pe - 9p 


P®ii + U, 


2 \ 3 


= vO 


p&u+pU 2 \ v 


;{pe ~ dp). 
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In view of (13.41) and (13.51) . we see that 

,1 — v ( p(Q ii + 022 + © 33 ) 1 1 — v ( pg\Ug\ 2 — 9p\U\ 2 


Tgu = d- 


Pe 




pe - dp 
Pe 


vd 


pQa + pU 2 


Pe 


(3.6) 


= 0 


1 “ V f p(Bll + 022 + © 33 ) 
3 1 Pe 


, 1 -^ \ pe\U9?-0p\U? 


Pe 


- vd 
Pe ~ dp 


Pe 


— —T 

— /ll 

Pe 


1 - v f pg\Ug\ 2 - dp\U \ 5 


Pe 


+ 


Pe 

pe - dp 
Pe 

pQa + pU 2 N 1 
Pe 


pe - dp 
Pe 


Therefore, applying Lemma 13.11 Lemma 13.21 and the estimate (1) and (2) above, we find 
that 


Tea < 


d{\+c„^m? + i_v(x + vw>)m , i + 


1 -Vm 


1 -VW) 


1 -Vm 


< 


1 + C9^/£(J) + Cy/£(t) 

1 -VW) 


This leads to 

(3.7) 


1 — v^W 


Lower bound estimate for Tgu can be derived analogously as 

(3.8) Tgu - 1 > -Cy/£(t). 

Combining (13.71) and (13.81) . we obtain 

\Teu - 1| < CVS (i). 

The case for i 7 ^ j is simpler. We first observe from the definition of Ggij that 


/ pgQgij + pUgiUgj 

V 2 


= dv 


p© i'j T pUiUj 


Tgij = —diOij + UiUj) - UgiUgj. 
pe 


Hence we have 

(3.9) 

Then we can proceed similarly to obtain the desired result. 
Lemma 3.4. Let £(t) be sufficiently small. Then we have 

( 1 ) \d a p e {x,t)\ < \ZW), 

(2) \d a Ug(x,t)\<C ]al £(t), 

(3) \d a Te(x,t)\ <C W \£{t). 


□ 


for some positive constant C| a |. 

Proof. The proof is almost identical to Lemma 13.21 We omit the proof. □ 








































20 


SEOK-BAE YUN 


Lemma 3.5. Let£(t) be sufficiently small. Then determinant of the temperature tensor If 
satisfies the following estimates: 

(1) |d“{det70}|, \d a {det% 8 }\ < C^£{£), 

(2) | det |, | det 70s \ > ^ , 

for a positive constant C independent of v. 

Proof. We recall the explicit formula for det 70 derived in the proof of Lemma 2.2: 

det 711 = TwTtTTsz — T^Tn — T 31 T 22 — 7/^ 7/3, 

det 70 = 70 ii70227033 ~ T^fTeu — Tg 31 Tg 22 — Tg 12 Te 33 . 

Then (1) follow from the direct application of the estimates on the derivatives of the macro¬ 
scopic fields in the preceding lemmas. To prove (2), we recall from Lemma 13.21 and Lemma 
13.31 that 


Tii = 1 + o(£(t)) ( i = 1,2,3), Tij = o (£(t )) (i ^ j). 

which leads to 


det T v , det Ten 


> 

> 


{i +om)} 3 

1 - o(£{t)) 

1 

2 


1 + 3{o(£(f))} 2 {! + o(£(f))} 


for sufficiently small £{t). 


□ 


3.2. Uniform estimate on the temperature tensor. Recall that the nonlinear pertur¬ 
bation T(/) contains inverse of the temperature tensor Tjf 1 : 


Q 


M 

ij 


1 Ptf{p,v-U e ,U e ,T e -\v) ( 1 

-ffiP -^ (w ,d e ,r»)- ex H'U 


Ue) T r e -\v 



Now, since 70 (and 70) contains v, rough estimates of its inverse may involve factors inversely 
proportional to v in it, which make it impossible to derive estimates uniform around v = 0. 
This is a serious problem considering that the v = 0 corresponds to the classical BGK model. 
In what follows, we will carefully investigate the temperature tensor 70 and show that the 
seemingly problematic \/v factor actually does not cause any harm. The key observation is 
that 70 is essentially equivalent to the temperature T under our assumptions on v. 


Proposition 3.1. Let — 1/2 < v < 1 . Define constant CO 1 and CO 2 by 
COi = min{l — v, 1 + 2v}, CO 2 = max{l — u, 1 + 2u}. 

Then the temperature tensor is comparable to the temperature in the following sense: 

CviT(x,t)Id < < Cv2T(x,t)Id. 

Furthermore, if £(f(t)) be sufficiently small, then 70 is invertible and 

Cv2 -Id < T-\x,t) < 3^-rId. 


T(x, t) 


T{x,t) 
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Proof. (1) We first observe from the definition of 7/ that 

(1 — v)pT + i/pQ n vpQ 12 vpQiz 

vpQ 21 (1 — v)T + upQ 22 vp<d 23 

I/0 3 1 vpQ 32 (1 - i/)T + i/p0 33 

= (1 — v)pTId + vpQ 


p7/ = 


■ / F(a:, u, f)|u — f7| 2 du/d + i/ j F(x,v,t)(v ~ U) <g> (v — U)dv. 
Jr 3 Jr 3 


Then a direct computation using 


fc T {(u -U)®(v- U)}k = {(« - U) ■ k} 2 


shows that for any k in R 3 

k T {pT v }k = -—-—- j f F(x,v,t)\v — U\ 2 dvX \k\ 2 + v f F(x, v, f){ (v — U) ■ k} 2 dv. 

3 L7r3 J Jr 3 

We split the estimate into the following two cases. When 0 < v < 1, we have 

k T {pT v }k > —— \k\ 2 f F(x, v, t)\v — U\ 2 dv. 

3 Jr 3 

In the case — ^ < v < 0, we apply Cauchy-Schwartz inequality to the second term to get 

k T {p%}k > ^ o ^ { f F(x,v,t)\v — U\ 2 dv\ \k\ 2 + is f F(x, v, t)\v — U\ 2 \k\ 2 dv 

\ Jr 3 J Jr 3 


3 

(1 + 2u) 


|fc | 2 f F(x, v, t)\v — U\ 2 dv. 
Jr 3 


Therefore, we have 

(3.10) k T {pT v }k > — min{l — u, 1 + 2z/}|fc | 2 j F(x, v, t)\v — U\ 2 dv, 

3 J R3 

or equivalently, 

(3.11) k T %k > min{l — z/, 1 + 2z/}|fc| 2 T, 

We then apply Lemma 13.II to compute 

T(x,t) = — [ F(x,v,t)\v — U\ 2 dv 
P Jr 3 

F(x, v , t)\v\ 2 dv — p\U\ 2 


1 


(3.12) 


1 - e(t) 


[ {p + ^JLf}\v\ 2 v~C£{t) 
Jr 3 


- i —f d\v\ 2 dv-\\f\\ L ~ f y/p\v\ 2 dv — C£(t) 


> 


1 - £(t) 


>3 -Cy/Itf) 

for some generic constant C. From (13.111) and (13.121) . we conclude that for any fixed —1/2 < 


v < 1 and for sufficiently small £(t), 7 / is invertible and 

1 


T~ < 


min{l — v, 1 + 2 i/} 


T~ l Id. 
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The proof for the upper bound is similar. □ 

Lemma 3.6. Let —1/2 < v < 1. Suppose £(f(t)) be sufficiently small. Then there exists a 
positive constant C v < oo such that 

X T {T- 1 }Y<C I/ {\\Xf + \\Y\\ 2 }, 

for X, Y in ffi 3 . 

Proof. By Proposition 13.II If is invertible under the assumption of the lemma. Moreover, 
Since 7f is symmetric, 7/T 1 also is symmetric. Therefore, we can compute 

|X t T" 1 T| = i \(X + Y) t T- 1 (X + Y)~X t T- 1 X-Y t T~ 1 Y\ 

< i \(X + Y) t T~\X + Y)\ + 1 \X T r~ l X\ + i \Y r T- l Y\ 

£ ™{T^rTM !l,x||2+l|y||2) - 

for any two vectors X and Y in R 3 . □ 

Similar result holds for Tg: 

Lemma 3.7. Let —1/2 < v < 1. Suppose £(f(t)) is sufficiently small. Then To is invertible, 
and there exists a positive constant C v < oo such that 

X T {T e - 1 }Y<C v {\\X\\ 2 + \\Y\\ 2 }. 

Proof. In view of (13.61) and (13.91) . we can write pTf as 

(pe\Ue\ 2 -ep\U\ 2 ) + pe-e P ^Id 

+ v8(p U ®U - pD) - v(p g Ug ®Ug - p s Dg), 

so that 

k T { Pe Te}k = dk{pT}k + {p e -e P )\k\ 2 + ^^{po\Ug\ 2 -Op\U\ 2 )\k\ 2 

+ vO {p(U ■ k) 2 - pk T Dk} - v {pg(Ug ■ k) 2 - pgk T Dgk} , 

for k £ M 3 . D and Dg denote the diagonal matrix with diagonal elements Uf,Uf,U 2 and 
u en u e 2 ’ U e 3 respectively: 

( U 2 0 0 \ f U 2 e 0 0 \ 

D = I 0 U 2 0 , Dg = 0 U 2 e 0 . 

V 0 0 U 2 ) V 0 o u 2 e ) 

Then, employing Lemma 4.1 and 4.3, we obtain 

k T {p g Tg}k > i min{l + 2u, 1 - v}\k\ 2 + (1 - 0)\k\ 2 + 0(£{t))\k\ 2 

> — min{l + 2u, 1 — r/}|fc| 2 , k £ M 3 , 

3 

for sufficiently small £(t). By virtue of Lemma l3.3l (1) 

k T {Tg}k > — min{l + 2v, 1 — v}\k\ 2 , 

The rest of the proof is similar to the proof of Lemma 13.61 □ 


PeTg = dpT + 


1-v 
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Lemma 3.8. Let — 1/2 < v < 1 . Suppose £(f(t)) is sufficiently small. Then there exists a 
positive constant C v ^ a < oo such that 

(1) X r {d a T~ 1 }Y < C v , a { ||X|| 2 + ||T|| 2 }, 

(2) X r {d a T~ g 1 }Y < a, Q {||X || 2 + ||F|| 2 }, 

for X, Y in R 3 . 

Proof. We have proved in Lemma T3.61 that 7/ is strictly positive definite for —1/2 < i/ < 1 
when £(t) is sufficiently small. Therefore, 7/ is invertible. Now, applying <9 to 7/7/T 1 = /, 
we see that <97/ {7/ -1 } + 7Z<9 {7/ -1 } = 0, and thus, 

(3.13) d{T- 1 } = T- 1 {d%}T~ 1 . 

Then the case |a| = 1 follows directly from this identity and Lemma 13.61 and Lemma |3.4I 
For general case, we recall 

(3.14) d°T u = ^ d^T.d 1 {T- 1 } 

I8I+M=M 

and use the induction argument. The proof for Te is almost identical. We omit it. □ 


3.3. Local existence. We first estimate the nonlinear term T(/). Note that, in contrast to 
the Boltzmann equation, we need to use the estimates on the macroscopic fields established 
in the previous section to control T(/) in the energy norm. 


Lemma 3.9. The bilinear perturbation T satisfies the following estimates: 


(1) 


( 2 ) 


9?W)gdv 


<C J2 \\ dai fhlJd a *f\\Ll\\h\\ Ll 


|ai| + |a 2 | 

<M 


+C J2 H^ 1 /II^JI^ 2 2 /llL ? l|5“ 3 /IU ? l|/ i ||L S , 


l“i| + l“ 2 |<|a| 
I< 83 I<I< 8 I 


+C J2 \\ dai f\\LlJ\d a2 f\\ L 2 v \\d a3 f\\ L i\\h\\ L 2, 


r i ,2 (f,g)fdv 
^3 (f,g,h)fdv 


|c*l| + |«2|-H“3l 

<M 

ri,2 (gj)fdv < CsupII^IIliJI/II 2 ^. 

3 {g,h,f)fdv 


r 3 (ff,/, h)fdv 

< Csup \\ghi sup ||/i||l2||/|| 2 2 iU . 

(3) ||ri, 2 (/,ff)/i + ri i 2 (5,/i)/i|| i 2 < C'sup |/i| sup ||/||z .2 ||^||i,2 

® v XV X 

l|r 3 (/, g, h)r + r 3 (g, f, h)r + T 3 (g, h, f)r\\ 

< Csup \h\ sup ||/|| L 2 sup \\g\\ Ll \\h\\ Ll . 


Proof. Recall that the T consists of Ti, T 2 and T 3 . We prove this lemma only for T 2 , because 
the proof for the remaining parts are similar. Utilizing macroscopic estimates established 
in the previous section, we find that there exists a polynomial P a .0, which is generically 
defined, such that 
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= C a ,p \P a ,p{v, dpe 1 dUg, dT v g)\ exp - UgYT^iy - Ug)^j 


< C a ,pP(v)ex p ^-{1 -o(£(f(t)))}'-Y+o(£(f(t))) 

< C a ,p,e exp ^-{l — o(£T(t))} Q- £ ) H 2 + o(£(f(t)))^J , 

where d denotes any of <9™ such that to < |a| and n < |/3|. Therefore, there exists a positive 
number a depending on a, f 3 and v such that 

1 


(3.15) 




dpQ M | < exp (-«H 2 ) 


(1) By (13.151) and Holder inequality, we see 

f \d%r 2 (f)g\dv < 

J R 3 


< 


< 


< 


< 


E 

H“al 
=M 

c E 

|c*i | + |ct21 + |Q 3 1 
= |a 3 | 

E 


|ai| + |“2| + |“3l ' 

= M 


/ R 3 


|ai| + |Q2| + |Q3l 

=M 


dpQ M {d ai /, 6i) (9 “ 2 /, | 

exp (~a\v\ 2 ) (d ai f, ei)(d a2 f,ej)gdv 
exp (-|^| 2 ) gdv) || 9 “ 1 /|| i 2 || 9 “ 2 /|| 


Ll 


'R 3 


c E ||exp(-«|r;| 2 )|| i2 J| ff || i2 „||9 Ql /IU ? ||9 a2 /IU ? 

|c*l |+|d£ 2 1 + |«*31 

= M 

C E \\d ai f\\Li\\d a 2 f\\ L z v \\g\\ L i. 


|ai| + |a2| + |a 3 | 

= M 


( 2 ) can be estimated similarly as 

/ r 2 { f , g)fdxdv < C [ \\ f \\ Li \\ 9 \\Li ( [ exp (~ a \ v \ 2 ) fdv\dx 
J Jr 3 \J r 3 J 


< C ||/|U;||5||L;||/||^d® 

Jr 3 

< C sup \\g\\ L 2 1|/|| 2 j 


(3) For $ £ L 2 , we have 


c [ 

ll/ll 

LlhWLlW'&WLldx 

Jr 3 


/* 

Csup 

M , 

/ ll/MMM*lh 

x,v 

j 

R 3 

Csup 

x,v 

M ( 

7 MllMlidx 

\Jr 3 . 


< C sup |r| sup \\f\\Ll\\g\\LzJ\®\\L% iV - 


X,V X 


Therefore, the duality argument gives 

l|r 2 (/, 5 )r || L 2 < C sup |r| sup ||/||l 3 IMIlJ ■ 


X.V X 
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□ 

From the estimates in Lemma 13.91 the following local existence theorem can be proved 
by standard arguments (See, e.g [Till PTTj l. 

Theorem 3.10. Let v be a fixed constant such that —1/2 < v < 1. Let Fq = go + y/fifo > 0 
and fo satisfies the conservation laws \2.5\) . Then there exists Mq > 0, T» > 0, such that if 
T* < and £(fo) < ^, there is a unique solution f(x,v,t) to the ES-BGK model \ 2. !J\) 
such that 

(1) The high order energy £(/(i)) is continuous in [0,T*) and uniformly bounded: 

sup £(f(t))<M 0 . 

o <t<T, 

(2) The distribution function remains positive in [0,T*): 

F(x , v,t) = p + y/jif(x, v, t) > 0. 

(3) The conservation laws \2.5\) hold for all [0,T*]. 

Proof. We consider the following iteration scheme. 

(3.16) d t F n+1 + v ■ S7 x F n+1 = {M v {F n ) - F n+1 } , 

where M.(F n ) is defined by 

M ^ n) - vmm exp G" - vn) ( ” - un) ) ■ 

p n , U n and Tf denote the local density, bulk velocity and the temperature tensor associated 
with F n = p + yffif n ■ With estimates on the nonlinear perturbation in Lemma [3.91 it is 
standard to prove the local existence (See [19] [41]). The only thing to be careful about is 
whether the temperature tensor Tff remains strictly positive definite for each n, so that the 
iteration scheme is well-defined in each step. But this follows directly from Proposition 13.II 
and Lemma 13.61 - 13.81 □ 


4. Global Existence 

Now, having all the necessary estimates at hand, the global existence can be established 
using standard arguments (See [13101]). We sketch the proof in this section. First, we need 
to recover the degeneracy of the linearized relaxation operator to obtain the full coercivity. 
For this, we define 

a(x,t)= / fy/fidv, bi(x,t)= / fvi^ffidv (i = 1,2,3), c(x,t) = / f\v\ 2 ^JIdv. 

J R 3 J R 3 Jr 3 

We also define a macroscopic projection P as follows: 

Pf = a(x, t)y/JI + ^2 b i ( x > t) v iVh + c(x, t)\v\ 2 ^JI. 

i 

Note that P is not identical to Pq but equivalent. Since L u {Pf} = 0 for — 1/2 < v < 1 by 
Corollary 12.11 we can split the linearized ES-BGK model (12.31) into the macroscopic part 
and the microscopic part as follows: 

{dt + v ■ V x }{Pf} = ~{d t + v • V,}{(/ - P)f} + L{(I - P)f} + F(/). 

We then expand the l.h.s and r.h.s with respect to the following basis (1 <i,j< 3): 

(4-1) {y/p,Viy/p,ViVjy/jL,v 2 ^Ji,Vi\v\ 2 ^Ji}, 
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and compare coefficients on both sides to obtain the following micro-macro system [19] : 


(4.2) 


dtd — £a T ha-i 


Ofbi + dxi d 
&xi bj + & X j bi 
&Xi bi + &tC 
&Xi d 


£abi + habii 

li 


ij i 

^ci H - K 


(* ^ jO 


hbcii 


for i,j = 1, 2, 3. Then, by carefully studying this system, we find that the macroscopic part 
can be controlled by the macroscopic part as follows (See [19]): 

(4.3) Y, {\\ d a a\\ L l + \\d a b\\ L l + \\d a c\\ L iJ<C £ ||d“(4 + MIUl- 

H<iv H<jv-i 

We slightly abused the notation on the r.h.s for the simplicity of presentation. On the other 
hand, we can bound i v and h u by the energy norm of / as 

£ l|3 Q (r + /oi \ L %<c v £ \\{i-p)d a f\\ Llv + c v y/W 0 £ \\d a f\\ Llv . 

|ce|<JV—1 |a|<JV |a|<A 

Combining this with (14.31) . we see that 

£ \\d a Pf\\l lv < i £ {\\d a a\\l lv + \\d a b\\l lv + 11^111^} 

|a|<iV ’ |a|<JV 

< c £ \\d a (I - P)f \\ 2 L 2 +CVm 0 £ ||5“/||| 2 , 

• ^ x,v • x ,v 

|a|<JV |a|<JV 

which implies 

(4.4) £ \\Pd a f\\l lv <c £ \\{I - P)d a f\\l ly 

|a|<iV ’ |a|<JV 

Therefore, Proposition 12.81 together with (14.411 and the equivalence of Pq and P imply the 
coercivity estimate for L v : There exists b v = S(u) >0 such that 

(4-5) £ (L„d a f,d a f) < -S v £ Wf{t)\\Li' V , 

|ct|<JV |ct|<JV 

when / is sufficiently small in the energy norm. We are now ready to derive the nonlinear 
energy estimates which enables us to extend the local solution into the global one. Let / be 
the smooth local in time solution constructed in Theorem 13.101 Taking <9“ on both sides of 
m, we obtain 

d t d a f + v ■ v,.a Q / = Ld a f + d a v(f). 


We then take inner product with d a f 

J t \\d a f\\h, v < (Ld a f,d a f) Llv + (d a T(f),d a f) Llv , 

and apply the coercivity estimates (14.51) together with the nonlinear estimates in Lemma 
13.91 to derive 

£ ii ^/ ii !^ < cvsimnm), 

|q|<A 


E° 
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where T>(f(t)) denotes 

*>(/(*))= E \wmi lv - 

M + |/3|<N 

We now turn to the general case involving the derivatives in the velocity variables. Applying 
dp to (12.31) , we get 

{d t + v ■ W x + v 0 }d a pf = d Pl v ■ S7 x d a p_pJ + dpPd a f + d a pY{f, /). 

We multiply dpf , integrate over x and apply Holder inequality with Lemma T3.9I to 


see 


E a a 


1 d_ 

2 dt 


l lv + voWWWh.* ^ C E 11^/IUs.JI^/IUs.. 


+ c\\d a f\\ Llv Wpf\\ Llv + cy/sjmnm)- 

Then, we split the first two terms in the r.h.s using Young’s inequality and gather relevant 
terms together to obtain 

+ C£(f{t))V{f{t)). 

Now, we observe that the r.h.s of Xq£|= m +i Ep can be controlled by the good terms of 

Cm E E^+CmY J E a 

|/31 =m a 

if Cm is sufficiently large. By good terms, we mean the production terms on the l.h.s. 
Therefore, we can find constants C m and S m inductively such that 


E g 


|a| + l/3|<JV, 
|/3| <m 


dt 


ot f || 2 | r || f\\2 

L 2 -h 0 m 


h\< c Ny /s(m)v(m). 


From this energy estimate, the existence of global solutions follows from the standard conti¬ 
nuity argument. Remaining part of the Theorem 1.1 can be established in the exactly same 
manner as in the classical BGK case [IT] . This completes the proof. 
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